On the phase transition of light in cavity QED lattices 
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Systems of strongly interacting atoms and photons, that can be realized wiring up individual 
cavity QED systems into lattices, are perceived as a new platform for quantum simulation. While 
sharing important properties with other systems of interacting quantum particles, here we argue 
that the nature of light-matter interaction gives rise to unique features with no analogs in condensed 
matter or atomic physics setups. By discussing the physics of a lattice model of delocalized photons 
coupled locally with two-level systems through the elementary light-matter interaction described by 
the Rabi model, we argue that the inclusion of counter rotating terms, so far neglected, is crucial 
to stabilize finite-density quantum phases of correlated photons out of the vacuum, with no need 
for an artificially engineered chemical potential. We show that the competition between photon 
derealization and Rabi non-linearity drives the system across a novel Z2 parity symmetry-breaking 
quantum criticality between two gapped phases which shares similarities with the Dicke transition 
of quantum optics and the Ising critical point of quantum magnetism. We discuss the phase diagram 
as well as the low-energy excitation spectrum and present analytic estimates for critical quantities. 

PACS numbers: 42.50 -o, 42.50 Pq, 73.43.Nq, 05.30.Rt 
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Introduction - Interaction between light and matter is 
one of the most basic processes in nature and represents 
a cornerstone in our understanding of a broad range of 
physical phenomena. In the study of strongly correlated 
systems and collective phenomena, light has traditionally 
assumed the role of a spectroscopic probe. The increas- 
ing level of control over light-matter interactions with 
atomic and solid-state systems 13] has brought forth 
a new class of quantum many-body systems realized on 
photon lattices [44la| where light and matter play equally 
important roles in emergent phenomena. 

The basic building block of such systems is the ele- 
mentary Cavity QED (CQED) system formed by a two- 
level system (TLS) interacting with a single mode of 
an electromagnetic resonator. When CQED systems are 
coupled to form a lattice, the interplay between photon 
blockade 17- T^] and inter-cavity photon tunnelling leads 
to phenomenology akin to those of Hubbard models of 
massive bosons as realized e.g. by ultracold atoms in op- 
tical lattices [2^] ■ The possibility of quantum phase tran- 
sitions of light between Mott-like insulating and super- 
fluid phases has stimulated a great deal of discussion re- 
cently 0-0; U " JJ I ■ The excitement about these systems 
stems from their potential as dissipative quantum simu- 
lators that provide full access to individual sites through 
continuous weak measurements [3]. 

While sharing important features with conventional 
condensed matter or atomic physics setups, systems of 
strongly correlated photons have their own unique prop- 
erties that ultimately derive from the nature of the fun- 
damental light-matter interaction. As photons can dis- 
appear by interacting with the matter field, their number 
is not conserved but rather fixed by the condition of ther- 



mal equilibrium. To describe this situation for a photon 
gas in equilibrium with either photonic or dipolar bath - 
such as in a blackbody - one says that photons have zero 



chemical potential [211 ] . From the point of view of bosonic 
Hubbard models this has rather dramatic consequences, 
as one would then require an external non-equilibrium 
drive in order to engineer non-trivial quantum many 
body states other than the vacuum [22|. For a Lattice 
CQED system however, as we will show in this Letter, 
this is remarkably not so. The light-matter interaction 
strength can play the role of an effective chemical poten- 
tial to stabilize finite density quantum phases of corre- 
lated photons out of vacuum. We find that the so-called 
counterclockwise terms in the elementary CQED Hamil- 
tonian, the Rabi model, typically neglected for quantum 
optical systems are relevant perturbations that explic- 
itly break the conservation of total number of excita- 
tions. This changes completely the nature of quantum 
criticality of lattice CQED systems from those based on 
Jaynes-Cummings model [4l4l5l]. leading to a Z 2 parity- 
breaking quantum phase transition where the two level 
systems polarize to generate a fcrroclcctrically ordered 
state and the photon coherence acquires a non-vanishing 
expectation value due to hopping. This novel quantum 
criticality, described by a delocalized super-radiant quan- 
tum critical point, shares some similarity with the Dicke 
phase transition of quantum optics and turns out to be 
in the universality class of the Ising model. 

Single Resonator - The elementary light-matter in- 
teraction between a photonic mode of a resonator and a 
TLS is described by the Rabi model [l[ 
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where uj r is the frequency of the resonator, uj q the qubit 
transition frequency, g the light-matter coupling strength 



and x 



In addition, depending on the spe- 



cific context, an extra term should be added to Eq. (fTJ) 
where the field appears quadratically, = D (a + a)) 2 . 
We will discuss its implications at the end of the Letter 
and first address the physics of the Rabi model (QJ and 
its lattice extension. When the coupling g is sufficiently 
smaller than the frequencies uj r , w q one can safely neglect 
processes where simultaneously atomic and photonic ex- 
citations are created, described by the counter-rotating 
terms <j + + a~ a. In this so called rotating-wave ap- 
proximation the Hamiltonian (Q]) reduces to the Jayncs- 
Cumming (JC) model, Hjc = w r a + uj q a + a~ + 
g (a' a" + a + a) used widely in discussions of CQED 
physics. While appropriate in many relevant cases, re- 
cent implementations of circuit QED [23j, |24| achieved 



coupling strengths g where the counter-rotating terms be- 
gin to show significant deviations from the expectations 
of the JC model [J [H-El . This is the so called "ultra- 
strong coupling" regime of parameters where g becomes 
a considerable fraction of uj r . 

Although the physics of the Rabi model has been 
widely studied and well- understood 32- 1H it is attract- 
ing renewed attention recently 35, 36l|. Here we will ex- 
amine the Rabi-Hubbard model as realized e.g. by a 
lattice of circuit QED cavities where each node is de- 
scribed by Hr. As we discuss below, this system forms 
a viable platform for studying non-trivial strongly corre- 
lated phases of light. Experimental efforts to fabricate 
on-chip photonic lattices of circuit QED systems are cur- 
rently underway 37j ■ Before we introduce the lattice, we 
consider the generalized Rabi model 

H g R[a, a^} = uj r a) a + uj q a + a~ + g (a^ o~ + a + a) + 

+g' (a t cr+ + a~a)(2) 

We would like to stress that we introduce this model 
to explore the role of counter-rotating terms in a con- 
trolled fashion. This model interpolates between the JC 
Hamiltonian for g' — and the standard Rabi Hamil- 
tonian for g' = g. In the following we will restrict our- 
selves to the resonant case lo t = lo q = luq. For g' = 
i.e. in the JC limit, the above model conserves the total 
number of excitations, Af — a + a + tj~ . The result- 
ing continuous U{1) symmetry allows an exact analytic 
solution of H g R in terms of dressed states of photons 
and TLS excitations, the polaritons. The ground-state 
shows an interesting evolution upon increasing the cou- 
pling g/ujQ, with an infinite series of level crossings for 
9c(n) = (V n + 1 + y/n) where the number of excita- 
tions increases from n to n+1, resulting in a characteristic 
staircase structure (see Fig. Q]). 

For an arbitrarily small g' , counter-rotating terms 
break the continuous U(l) symmetry down to a discrete 
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Z2 group associated with parity V = e' 
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FIG. 1: Evolution of the single site generalized Rabi Model 
ground state properties upon increasing the strength of 
counter-rotating terms g'/g- Top panels: number of exci- 
tations and parity of the ground state as a function of g/uio 
and for different g'/g. Bottom panels: (left) decay of level 
splitting A in the ultra-strong coupling regime and (right) 
probability of having n polaritons in the Rabi ground state 
(inset) and its evolution with g'/g = 0.0, 0.25, 0.5, 0.75, 1.0. 



unitary operator the photon field a and the TLS op- 
erator a x transform respectively as V^aV = —a and 
'P^o'x'P = —cr x , from which the invariance immediately 
follows, namely [P,W g R] = 0. A direct consequence of 
the parity symmetry is that, while (a) = (<r x ) = in 
the ground state of HgR much as in the JC limit, the 
photon field in the Rabi ground state is squeezed, i.e. 
(a 2n ) 7^ 0. While the discrete Z2 symmetry prevent a 
full closed-form solution, the model ([2]) in the Rabi limit 
g' = g has recently been shown to be nevertheless in- 
tegrable (35[. Important features of this exact solution 
that will be relevant for our discussion below are that (i) 
no level crossing between states of different parities can 
occur as a function of g/co (note that g' = g), which in 
turn implies that the ground state of the Rabi model re- 
mains an even parity state for any g/wo, (ii) the ground 
state and the first excited state are quasi-degenerate in 
the ultra-strong coupling regime g > ojq ■ 

To get further insight into the structure of the Rabi 
ground state we numerically diagonalize the Hamilto- 
nian ([2]). In Fig. [T] (top panels) we plot the number of ex- 
citations and the parity of the ground state as a function 
of g/ujQ for different values of g'/g. Upon increasing the 
strength of counter-rotating terms the JC plateaux are 
gradually smeared out. Though the parity remains well- 
defined, the evolution with g' reveals multiple crossings 
between eigenstates switching the parity of the ground 
state, ultimately resulting in an even parity ground state 
when g' = g. We also plot (bottom panel, Fig. [T]) the 
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scaled level splitting A/uj q between the ground state and 
the first excited state, which vanishes as A ~ g-Hs/^o) 
for large g/ujQ in agreement with degenerate perturba- 
tion theory. In contrast, the gap to the next energy level 
stays of order one (not shown) at large g/uj . Polari- 
tonic dressed states are not anymore exact eigenstates 
for g' ^ 0. Turning on g' results in a broadening of the 
eigenstates of the generalized Rabi model when projected 
on the polaritonic eigenstates of the JC model, its peak 
shifting as a function of g' resulting in a ground state 
containing a finite number of excitations (right panel of 
Fig. P . 

Lattice model of interacting atoms and photons - We 
now come to the main subject of this Letter, which is 
the physics of the Rabi- Hubbard model, a model of itin- 
erant photons hopping between neighboring resonators 
and interacting on-site with a TLS according to the local 
Hamiltonian ©. The full many-body Hamiltonian for 
this system reads 



H = — J ^2 a R a R ' 

<RR') 



R 



Hg R [au,o-^}. (3) 



We stress here that, with an eye on possible future exper- 
iments on a circuit QED architecture, we do not include 
any chemical potential to tune the density of excitations 
in the ground state. The goal here is to exploit the spon- 
taneous polarization of the Rabi vacuum that emerges 
when cither the light-matter interaction strength g or 
the hopping strength J is increased. 

We start considering first the g' = JC limit and study 
the phase diagram in the g — J plane [42j • In the absence 
of hopping, the ground state is an exact dressed state of 
polaritons. A gapped and incompressible Mott Insulating 
(MI) phase of polaritons survives at finite hopping until 
a critical value of J is reached. The phase boundary 
J c (g) features characteristic Mott lobes (see Fig. [3]), a 
legacy of the level crossings of the single site JC model 
discussed above. For hopping strengths larger than J c 
the system is in a superfluid (SF) compressible phase with 
gapless excitations associated to phase fluctuations of the 
U(l) order parameter. It is now well-established that the 
JC lattice model is in the same universality class as the 
Bose-Hubbard model @, EHiS, [1M [3^] . Crucially for our 
purpose here, the experimental realization of this MI-SF 
quantum phase transition requires an external driving 
or a suitably engineered chemical potential in order to 
counter-balance photon losses into the vacuum. 

We now argue that the inclusion of counter-rotating 
terms in Eq. © has a dramatic effect on the above 
physics. The roots of this can be traced back to the 
single resonator limit. As discussed above, the counter- 
rotating terms leave the system with a discrete Z2 sym- 
metry associated to parity. Photon hopping in ((3]) can 
trigger a spontaneous breaking of this parity symmetry 
above some critical coupling J c {g), toward a phase where 



both ( or) 7^ and ( er R ) ^ 0. As the broken symmetry 
is discrete, this quantum phase transition is fundamen- 
tally different from the JC one. Indeed it can be seen as 
a delocalized super-radiant quantum critical point rem- 
iniscent of the multi-mode Dicke transition of quantum 
optics. In order to see that a non-zero J favors ordering, 
we start from the full Hamiltonian © and notice that 
photons can be integrated out exactly in an imaginary- 
time action formalism to obtain an effective model for 
the TLSs only. 
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The result of this calculation 
veal that photon mediates an effective Ising-likc cou- 
pling between TLS which is retarded and long-range, 
J r-r<( t ) = -9 2 /2(r T x R (r).T R ,(0)). The scaling with 
g implies that at sufficiently large g/coo and for finite 
J, a ferromagnetically ordered Ising phase emerges with 
(<r R ) 7^ 0. Further insight into this emerging Z2 degree 
of freedom are obtained from the single site limit. As we 
discussed, at large g/uio the ground state and the first ex- 
cited state are almost degenerate, with an exponentially 
small splitting and a gap to the next level which stays of 
order one. These two states |±) have opposite parity and 
can be thought as eigenstates of an effective pseudospin 
1/2 degree of freedom, £ R . In addition we notice that 
the photon operator a R does not couple states with same 
parity. Its expression in the restricted |±) subspace reads 
a R — > (3 £ R + 7 £ R , where the dependence of the coeffi- 
cients /?, 7 on the coupling g can be obtained numerically. 
In the limit g/u>o 3> 1 one can analytically show that a 
linear scaling holds (3 = 7 ~ g/^o- Armed with these 
results we can rewrite the Rabi- Hubbard Hamiltonian as 



H 



eff 



£ J*£ R £ R ,+.F5] R £ R , + -£E R (4) 



(RR'j 



R 



where the couplings J x ' v = J (7 ± /3) 2 /2 depend on g as 
shown in Fig. [5J This effective model describes a pseu- 
dospin anisotropic XY model in a longitudinal magnetic 
field A/2, which is known to display a quantum phase 
transition toward a Z2 broken symmetry phase which is 
in the Ising universality class for any finite anisotropy, 



jx j V |4jJ] _ The effective psucdospin description high- 
lights once more the differences between the Rabi and 
the JC case. Indeed here both the disordered and the or- 
dered phases are gapped except right at the critical point 
where the gap is expected to vanish as a power-law. 

Mean Field Phase Diagram and Fluctuations - We 
now use Gutzwillcr mean field theory to confirm the gen- 
eral picture we have drawn for the transition. By de- 
coupling the hopping term in ([3]) we reduce the origi- 
nal lattice problem to an effective single site problem, 
Hef f bP] = Hi oc — if; (at + a) in a self-consistent field 
ip = Z J(a)ip. By expanding the energy to second order 
in ip we can get the mean field phase boundary J c (g), 
above which a parity symmetry broken phase emerges 
with both (o a ) ^ and (er R ) ^ 0. In Fig. [3J we 
plot the mean field phase boundary in the J, g plane for 



4 




FIG. 2: Top panels: (Left) Scaled ferromagnetic coupling 
J x ' y I J of the effective spin model as a function of g/u>o. While 
J x ~ (g/ujo) 2 at ultra-strong coupling, the coupling along y 
is vanishingly small. (Right) Mean-field order parameter for 
the effective spin model. Bottom panel: Band of spin-wave 
excitations above the ground-state. 



different values of g'/g from the JC to the Rabi limit. 
The "Mott lobes" for g' = are gradually suppressed 
as the ratio g'/g is increased. For intermediate val- 
ues a residual lobe structure remains, which reflects the 
level crossings already discussed in the single-site prob- 
lem. However we stress that no Mott insulator exists 
for any finite g'. Further insight on the transition can 
be gained from the effective spin Hamiltonian A 
linearized fluctuation analysis g^ives a critical coupling 
J c = A/ ((3 + 7) 2 ~ wjj e" 2 ^/^ ) /4a 2 which agrees with 
the numerical results in the large g/uJo regime (see fig- 
ure [3J. In addition the effective Hamiltonian also gives 
access to the spectrum of low-lying excitations Uk, plot- 
ted in Fig. [21 that as expected is gapped on both sides 
and vanishes in a power-law fashion at the transition, 
Eg ~ | J — Jcl 1 ' 2 . The low-energy spectrum is gapless at 
the critical point, with a linear dispersion cjk = c |k|. 

Discussion - The physical picture we have drawn from 
our analysis of the lattice Rabi model reveals a striking 
feature of hybrid systems made of atoms and photons. 
Due to the nature of the fundamental light-matter inter- 
action, which allows non-trivial vacuum fluctuations, no 
external driving forces or artificially engineered chemi- 
cal potentials are in principle required to stabilize finite 
density quantum phases of correlated atoms and photons. 
Rather it is the coupling between matter and light that 
will trigger this non-trivial vacuum polarization. An im- 
portant question for a possible experimental realization 
e.g. on a circuit QED platform concerns the stability of 
the above picture against photon leakage that is an in- 
herent feature of any quantum optical system. Physical 




J/co 



FIG. 3: Mean field phase diagram for the generalize Rabi 
lattice model for different values of g'/g. For the Rabi case, 
g' — g, we compare with the mean field phase diagram for the 
effective spin model (left bottom panel). In the inset, we show 
the decay of the critical coupling J c at ultra-strong coupling 



intuition would suggest that at least for a small coupling 
to a low-temperature photonic bath the ordered phase 
would be protected by the discrete nature of the Zi sym- 
metry. However an in-depth study of the phase diagram 
and a full understanding of quantum criticality in the 
open system limit is an important fundamental problem 
that we leave to future investigation. 

We now briefly discuss the role of a generic A 2 term 
H A 2 — D (a + o)) 2 in the picture that emerges from the 
above discussion. Recently it has been argued that such 
a term, generally assumed to be small, can become rel- 
evant in certain cavity QED realizations of the single- 
mode Dicke Model, where an ensemble of many TLSs 
is coupled to a single mode of a cavity. Indeed when 
the coupling D scales sufficiently fast with light mat- 
ter interaction g, D > g 2 /ui q , the super-radiant criti- 
cal point disappears. While this condition is realized in 
cavity QED setups with real atoms coupled via electric 
dipolc and results in so called no-go theorems, the sit- 
uation with some circuit QED implementations, where 
TLSs couple capacitively to the resonator, is currently 
subject of scientific debate [HI, Hi}. We note that in con- 
trast to the single-mode Dicke model, in our system the 
Z2 parity symmetry breaking emerges from a non-trivial 
interplay between derealization of photons through hop- 
ping and local light-matter interactions. As a result the 
critical boundary J c (g) can be accessed by increasing the 
hopping strength J, at fixed (and even moderate) light- 
matter coupling. While the inclusion of H ^42 in our lat- 
tice Hamiltonian may change quantitatively the shape 
of the phase boundary [43| . especially in the ultrastrong 
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coupling regime, it is rather an issue of the specific im- 
plementation that will determine the ideal architecture 
to realize the Rabi phase transition in an experimental 
system. Finally we note that circuit QED implementa- 
tions can be engineered where the A 2 term is irrelevant. 
This is the case, for example, of flux qubits inductively- 
coupled to resonators [I?} , a setup that in principle HH is 
ideally suited to access the ultra-strong coupling regime 
and that has been recently explored experimentally in 
single circuit QED units j23j . 

Conclusions - In this work we have explored the 
physics of itinerant photons hopping between neighbor- 
ing resonators of a lattice of CQED systems. We have 
studied its equilibrium phase diagram as a function of the 
atom-photon coupling g and shown that this system dis- 
plays a novel parity symmetry breaking quantum phase 
transition, belonging to the Zi Ising universality class, 
between two gapped phases. Simultaneously, the pho- 
tonic degrees of freedom acquire a non-vanishing expec- 
tation value, displaying a delocalizcd superradiant phase 
above a critical hopping. 

Note added - The following reference that is rele- 
vant to the discussion in this manuscript was brought to 
our attention after the publication of this work. 

Acknowledgements - We thank David Huse for valu- 
able and stimulating discussions. This work was sup- 
ported by the National Science Foundation through the 
Princeton Center for Complex Materials under grant no. 
DMR-0819860 and by the Swiss NSF through grant no. 
PP00P2-123519/1. 



[2 
[3; 

[4: 
[5; 

[T 

[8 

[9; 

[10 

in 



and L. C. L. 



S. Haroche and J. Raimond, Exploring the Quantum: 
Atoms, Cavities, and Photons (Oxford University Press, 
USA, 2006), 1st ed., ISBN 0198509146. 
Y. Yamamoto, F. Tassone, and H. Cao, Semiconductor 
Cavity Quantum Electrodynamics (Springer, 2000), 1st 
ed., ISBN 3540675205. 

R. J. Schoelkopf and S. M. Girvin, Nature 451, 664 
(2008). 

M. J. Hartmann, F. G. S. L. Brandao, and M. B. Plenio, 
Nat Phys 2, 849 (2006). 
A. D. Greentree, C. Tahan, J. H. Cole, 
Hollenberg, Nat Phys 2, 856 (2006). 
D. G. Angelakis, M. F. Santos, and S. Bose, Physical 
Review A 76, 031805 (2007). 

D. Rossini and R. Fazio, Physical Review Letters 99, 
186401 (2007). 

M. Aichhorn, M. Hohenadler, C. Tahan, and P. B. Lit- 
tlewood, Physical Review Letters 100, 216401 (2008). 
N. Na, S. Utsunomiya, L. Tian, and Y. Yamamoto, Phys- 
ical Review A 77, 031803 (2008). 

I. Carusotto, D. Gerace, H. E. Tiireci, S. De Liberato, 
C. Ciuti, and A. Imamoglu, Physical Review Letters 103, 
033601 (2009). 

S. Schmidt and G. Blatter, Physical Review Letters 103, 
086403 (2009). 



[12 

[13; 
[14; 

[is: 
[is; 
[17; 
[is; 

[19 

[20 
[21 

[22 
[23 

[24 

[25 

[26 
[27 

[28 

[29 

[30 

[31 

[32 

[33 

[34 

[35 
[36 
[37 
[38 

[39 

[40 

[41 

[42 



Fazio, D. Gerace, 
Imamoglu, Physical 

Physical Review A 

Nat Phys 8, 292 



J. Koch and K. Le Hur, Physical Review A 80, 023811 

(2009) . 

P. Pippan, H. G. Evertz, and M. Hohenadler, Physical 
Review A 80, 033612 (2009). 
A. Tomadin, V. Giovannetti, R. 
I. Carusotto, H. E. Tiireci, and A. 
Review A 81, 061801 (2010). 
R. O. Umucahlar and I. Carusotto, 
84, 043804 (2011). 

A. Houck, H. E. Tiireci, and J. Koch, 
(2012). 

L. Tian and H. J. Carmichael, Physical Review A 46, 
R6801 (1992). 

A. Imamoglu, H. Schmidt, G. Woods, and M. Deutsch, 
Physical Review Letters 79, 1467 (1997). 
S. Rebic, S. M. Tan, A. S. Parkins, and D. F. Walls, 
Journal of Optics B: Quantum and Semiclassical Optics 
1, 490 (1999). 

I. Bloch, J. Dalibard, and W. Zwerger, Reviews of Mod- 
ern Physics 80, 885 (2008). 

L. D. Landau and E. Lifshitz, Statistical Physics, Third 
Edition, Part 1: Volume 5 (Butterworth-Heinemann, 
1980), p. 171, 3rd ed., ISBN 0750633727. 
J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Nature 
468, 545 (2010). 

T. Niemczyk, F. Deppe, H. Huebl, E. P. Menzel, 
F. Hocke, M. J. Schwarz, J. J. Garcia-Ripoll, D. Zueco, 
T. Hummer, E. Solano, et al., Nat Phys 6, 772 (2010). 
P. Forn-Diaz, J. Lisenfeld, D. Marcos, J. J. Garca-Ripoll, 

E. Solano, C. J. P. M. Harmans, and J. E. Mooij, Physical 
Review Letters 105, 237001 (2010). 

V. Hussin and L. M. Nieto, Journal of Mathematical 

Physics 46, 122102 (2005). 

J. Larson, Physica Scripta 76, 146 (2007). 

J. Hausinger and M. Grifoni, New Journal of Physics 10, 

115015 (2008). 

T. Werlang, A. V. Dodonov, E. I. Duzzioni, and C. J. 
Villas-Bas, Physical Review A 78, 053805 (2008). 
J. Bourassa, J. M. Gambetta, A. A. Abdumalikov, 
O. Astafiev, Y. Nakamura, and A. Blais, Physical Re- 
view A 80, 032109 (2009). 

J. Casanova, G. Romero, I. Lizuain, J. J. Garca-Ripoll, 
and E. Solano, Physical Review Letters 105, 263603 

(2010) . 

F. Beaudoin, J. M. Gambetta, and A. Blais, Physical 
Review A 84, 043832 (2011). 

R. Graham and M. Hohnerbach, Zeitschrift fur Physik B 
Condensed Matter 57, 233 (1984). 

M. Kus and M. Lewenstein, Journal of Physics A: Math- 
ematical and General 19, 305 (1986). 
H. G. Reik and M. Doucha, Physical Review Letters 57, 
787 (1986). 

D. Braak, Physical Review Letters 107, 100401 (2011). 

E. Solano, Physics 4 (2011). 

A. A. Houck (Private communication). 

A. Tomadin and R. Fazio, Journal of the Optical Society 

of America B 27, A130 (2010). 

M. Hohenadler, M. Aichhorn, S. Schmidt, and L. Pollet, 
Physical Review A 84, 041608 (2011). 
S. Schmidt and G. Blatter, Physical Review Letters 104, 
216402 (2010). 

S. Sachdev, Quantum Phase Transitions (Cambridge 
University Press, 2011), 2nd ed., ISBN 0521514681. 
While it should now be clear that JC model is insufE- 



6 



cient to describe the ultra-strong coupling limit g' > ui q , 
its discussion will be helpful in identifying the new fea- 
tures that arise with respect to earlier studies that con- 
sider the JC phase diagram in the fi — J plane, where fj, 
is an artificially introduced chemical potential for polari- 
tons @, 1, @, Gl H, Hi . 

[43] See Supplementary Material for details on the effective 
spin model and the role of A 2 term. 

[44] We expect the Gutzwiller mean field results to be ac- 
curate in the regime J <C wo- The behaviour at large 
hoppings J ~ cjo requires a different approach as the 
tight-binding model (3) itself becomes questionable in 
this regime, where a significant mixing of single-resonator 
modes occur. 

[45] O. Viehmann, J. von Delft, and F. Marquardt, Phys. 

Rev. Lett. 107, 113602 (2011). 
[46] P. Nataf and C. Ciuti, Nat. Commun. 1, 72 (2010). 
[47] P. Nataf and C. Ciuti, Phys. Rev. Lett. 104, 023601 

(2010) . 

[48] M. Devoret, S. Girvin, and R. Schoelkopf, Annalen der 

Physik 16, 767 (2007), ISSN 1521-3889. 
[49] H. Zheng, Y. Takada, Physical Review A 84, 043819 

(2011) . 



Supplementary Material for "On the phase 
transition of light in cavity QED lattices" 

Integrating out the Photons: Effective Spin-Spin 
Interactions 



An interesting perspective on the physics of a model 
of localized spins and itinerant photons can be obtained 
by integrating out the photonic degrees of freedom. This 
can be done exactly as the photons always enter quadrat- 
ically. Let us first rewrite the quadratic part of the action 
associated to the Rabi model (we consider this case for 
simplicity, extension to the generalized Rabi is straight- 
forward). After diagonalizing the hopping in Fourier 
space we get 

So= dr a k (r) (d T - a; k ) a k (r) = (5) 
^ k 

/ dT^M^g^iT-rra^r') (6) 

J k 

where the photonic dispersion reads w k = Q — 
2J^ Q cos k a and we have introduced the bare photon 
Green's function 

G k (r) = -(T T a k (r)al(0)) (7) 

The coupling between the photon and the spin is linear 
in the field x, = a, + a,- 



S g = g V Xi{T)af(r) 

Jo 



(8) 



which after Fourier transforming reads 



= g f drJ2 Ok(r) t£ k (r) + 4(r) o£(Y) . (9) 
Jo k 

The integration over the photon field can be done exactly 
by a simple shift of variable and the result reads 



Z= / DaDaDae- {So+s <> +Su 



(10) 

with a spin-only effective action of the form 
Seff = S loc + g 2 [ drdr' £ <^ k (r) Q*(t - t') o£{V) 

3 k 

(11) 

where S\ oc = oj q ^ a t is proportional to the total inver- 
sion of the TLSs, which acts as a transverse field. We 
can rewrite this effective action as following 

S eff = S loc + f drdr 1 ^-kW Jk (r - r') a£(r')[12) 

J k 

where the effective spin-spin interaction is given in terms 
of the Green's function of the photon field 



(13) 



We notice this coincides with the Matsubara Green's 
function of the photon displacement 

Jk(r) = -^(T r x k (r).T_ k (0)) (14) 

a result which is also expected by looking at the order by 
order perturbation theory. 



The Rabi Lattice model with the A term 

In this section we discuss the modification of the phase 
boundary of the Rabi Lattice model when the A 2 term 
of the form 



R 



(15) 



dr J2 x^{t)o x _^{t) 
k 



is added. For the sake of simplicity we consider the case 
g' = g, hence the full Hamiltonian reads 

h=-j a R aRi +Y. UR ^ a ^ +HA2 ( 16 ) 

<RR') R 

where Hr is the hamiltonian of the single site Rabi 
model. We start noticing that, at a general level, this 
extra term does not break explicitly the Z2 symmetry 
of the problem, hence it doesn't a priori rule out the 
existence of the phase transition. At the same time we 
expect this term to modify the phase boundary in such 
a way to favour the disordered (symmetric) phase. 
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We notice at this point that, contrarily to the case of 
single mode Dicke Model, our lattice Hamiltonian cannot 
be solved exactly in general. In order to proceed it is use- 
ful to perform a canonical transformation on the photon 
fields to eliminate H^2 . This can be done by expressing 
the photon operators a R ,a R on each lattice site in terms 
of a new set of degrees of freedom 6 R , 6 r 



ihBb 



R 



By choosing 



or = coshO 6 R 



a R = sink® bfi + coshO 6 R 



tank 26 = -2D/ (ui r + 2D) 



(17) 
(18) 



we can eliminate the A 2 term and obtain a local Hamil- 
tonian which is again of the Rabi form, yet with renor- 
malized photon parameters 



H, 



rabi 



Ul r 6d b 

R 



R R 



^R+5cr R (6 R + 6 R ) (19) 



where 



r = ^/uJ 2 +4D 



U! r 



9 = 9 



We further notice that the hopping term between differ- 
ent resonators gets also renormalizcd and in addition a 
pair-hopping term arises so that the hopping Hamilto- 
nian reads in terms of the new fields 



H hopping — — J ^ (^r^R' + ^ c 

(RR'J 



(20) 



'J pair ^ (&R &R' + & R 6 R , ) (2'\) 
<RR') 



where the renormalized hopping read respectively 



J = J- 



UJ r + 2D 



J,, 



-J- 



2D 



(22) 




(0 /<D 




FIG. 4: Left Panel: Renormalized light-matter coupling g/Q r 
and TLS frequency uj q /Cj r as a function of the bare coupling 
g/uir for different values of the parameter a. We notice that 
for small g/uir the effects of renormalization are negligible for 
any a and become more relevant for larger values of g/cj r - 
Right Panel: Mean Field Phase diagram for different values 
of the parameter a. 



We start noticing that, quite generically for D > 0, the 
pair hopping term satisfies the condition 



pair | 



2D / UJ r 

1 + 4D/lj, 



J J 

< 

Ulr Ul. 



from which we conclude that if we assumed in the origi- 
nal problem J < ui r (an assumption which is behind the 
Hamiltonian [16]) then, within the same accuracy, we can 
drop the pair hopping terms in the transformed hamil- 
tonian. From the above discussion we conclude that the 
effect of the A 2 term can be fully encoded in a new lattice 
Rabi model with renormalized parameters. Whenever 
the coupling D can be made independent of the light- 
matter coupling g, a situation that can arise in certain 
circuit QED architectures, the effect of the A 2 term is 
just a trivial uniform shift of the phase boundary in the 



g vs J plane. The opposite situation is more tricky. We 
will assume in the following a coupling D that depends 
— < — (23) on the light-matter interaction q as D = a— and dis- 

cuss the physics of the model for different values of a. In 
order to gain insight, it is useful to look at how the cou- 
pling D renormalizes the photon parameters. In Figure 
H] we plot for different values of a the renormalized ef- 
fective light-matter coupling g/uj r and the renormalizcd 
effective TLS frequency ui q /Cj r as a function of g/uj r . We 
see that in general the role of a finite a is (i) to scale 
down the effective light-matter interaction particularly 
in the ultrastrong coupling regime g/ui r ^> 1, and (ii) 
to renormalize down the TLS frequency ui q /u) r <C 1, to- 
ward a regime of large detuning. We notice that the two 
effects compete with each other, as the former reduces 
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the tendency to order while the latter enhances it. This 
competition, in the case of single-mode Dicke model, has 
a rather drastic consequence for sufficiently large a - the 
disappearance of the supcrradiant critical point. How- 
ever, while the phase transition in the Dicke model is 
driven only by the local physics, namely g vs. oj q , in 
the Rabi model increasing the hopping strength J is an- 
other route to access the super-radiant phase. For this 
we find that the hopping strength is only weakly affected 
by the renormalization due to D. It is easy to see using 
Eq. and the definition of D that J/uj r ~ J/u) r , even 
at ultra-strong coupling. In conclusion, we expect the 



critical point to be accessible even in presence of the A 2 
term. While a thorough discussion of this issue will be 
presented elsewhere, here we give support to this claim by 
solving the Rabi Hubbard model with the A 2 term within 
the Gutzwiller mean field approximation As we see from 
the right panel of Figure [J] the shape of the phase bound- 
ary is only quantitatively affected by finite a, particularly 
in the ultra-strong coupling regime, but the overall pic- 
ture is not changed, namely a hopping driven quantum 
phase transition between a disordered normal phase and 
a dclocalizcd super-radiant phase where the TLSs order 
fcrromagnetically. 



